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We propose a feedback scheme to control the vibrational motion of a single trapped particle 
based on indirect measurements of its position. It results the possibility of a motional phase space 
uncertainty contraction, corresponding to cool the particle close to the motional ground state. 



I. INTRODUCTION 

In recent years there has been an increasing interest on trapping phenomena and related cooling techniques [Q. 
Some years ago it has been shown that, using resolved sideband cooling, a single ion can be trapped and cooled 
down near to its zero-point vibrational energy state |^] and recently, analogous results have been obtained for neutral 
atoms in optical lattices || . The possibility to control trapped particles, indeed, gave rise to new models in quantum 
computation Q, in which information is encoded in two internal electronic states of the ions and the two lowest 
Fock states of a vibrational collective mode are used to transfer and manipulate quantum information between them. 
It may happen however, that a trapped ion that is a favorable candidate for quantum information processing since 
it posseses a hyperfine structure with long coherence times (as for example 25 Mg + ||), is not suitable for resolved 
sideband cooling. In such cases it may be helpful to have an alternative cooling technique, which can be applied when 
resolved sideband cooling is impractical to use. 

In this paper we present a way to control the motion of a trapped particle, which is able to give a significant 
phase-space-localisation. The basic idea of the scheme is to realize an effective and continuous measurement of the 
position of the trapped particle and then apply a feedback loop able to decrease the position fluctuations. Due to the 
continuous nature of the measurement and to the effect of the trapping potential coupling the particle position with 
its momentum, feedback will realize an effective phase-space localisation. 

With this respect there are some analogies between the present method and resolved-sideband stimulated Raman 
cooling H , which can be viewed as a sort of feedback scheme. In fact, one of the two Raman lasers performs an effective 
measurement of the vibrational number by changing the particle internal state only if it is an excited vibrational state. 
The second Raman laser performs instead the feedback step, because it puts the particle back in the initial internal 
state, after having removed a vibrational quantum. The feedback scheme proposed here measures the particle position 
rather than its energy and tries to achieve cooling as phase-space localisation, using the particle oscillatory motion to 
mix position and momentum quadratures. 

A second analogy is given by the fact that the proposed method needs a Doppler pre-cooling stage, as it happens 
for resolved sideband cooling. In fact, the effective trapped particle position measurement is realized only in the 
Lamb-Dicke regime, i.e. when the recoil energy is much smaller than the energy of a vibrational quantum, which can 
be obtained only when the particle has undergone a preliminary cooling stage. Our scheme will provide therefore 
further phase space localisation and cooling. 

The paper is organized as follows. In section II we show how to realize the indirect continuous measurement of the 
position by coupling the trapped particle with a standing wave. In section III we shall introduce the feedback loop, 
in section IV we shall study the properties of the stationary state in the presence of feedback and section V is for 
concluding remarks. 
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II. CONTINUOUS POSITION MEASUREMENT 



We consider a generic particle trapped in an effective harmonic potential. For simplicity we shall consider the 
one-dimensional case, even if the method can be in principle generalized to the three-dimensional case. This particle 
can be an ion trapped by a linear rf-trap Q or a neutral atom in an optical trap [[U^]. Our scheme however does 
not depend on the specific trapping method employed and therefore we shall always refer from now on to a generic 
trapped "atom". 

The trapped atom of mass to, oscillating with frequency v along the x direction and with position operator x = 
xo(a + a^), xq — (fi^mv) 1 / 2 , is coupled to a standing wave with frequency cob, wave- vector k along x and annihilation 
operator b. The standing wave is quasi-resonant with the transition between two internal atomic levels |+) and |— ). 
The Hamiltonian of the system is || 

H = ~^~ a z + fiva) a + hujbb^b + ihe(a + + <J-){b — b^) sin (kx + 4>) , (1) 

where <j z = |+)(+| — |— )(— |, c± = |±)(=p|, and e is the coupling constant. In the interaction representation with 
respect to H = Tilu (b'b+ ^f), where ui ~ iOb will be specified later, and making the rotating wave approximation, 
this Hamiltonian becomes 

H = ^y°z + Tiva)a + KStfb + ihe(a + b - cr_6 t ) sin (kx + 4>) , (2) 

where A = loq — lu and 6 = ujb — u) are the atomic and field mode detuning, respectively. 

There are now two different ways for realizing an effective continuous measurement of the atom position and we 
shall describe them separately, even if they present many similarities. 



A. The Resonant case 



We consider the case when the standing wave is perfectly resonant with the |+) <-> |— ) transition, i.e., LUb = wq. It is 
therefore convenient to choose the frequency of the rotating frame lu = ujb = ujq in this case, so that both detunings are 
equal to zero. Moreover we shall consider the case of a very intense standing wave, so that it can be treated classically, 
that is, b can be replaced by the c-number j3. Choosing the phase of the field such that (3 = —i\/3\, Hamiltonian (||) 
becomes 

H = Tiva) a + he\(3\a x sin {kx + 4>) , (3) 

where a x = a + + cr_. If we finally set the spatial phase cf> = (i.e. the atom is trapped near a node of the classical 
standing wave) and assume the Lamb-Dicke regime, we can approximate the sine term at first order and get Jlfj| , pl| 

H = fiva)a + hx(J x X , (4) 

where x — 2e|/3|fcxo is the effective coupling constant between the internal and the vibrational degrees of freedom, 
and X — (a + a^)/2 is the dimensionless position operator of the trapped atom. This Hamiltonian shows how one can 
realize an effective measurement of the atomic position. In fact, the atom displacement away from the electric field 
node increases the probability of electronic excitation, and hence displacements can be monitored by means of the 
atomic fluorescence. Therefore, the two-level (sub)system can be used as a meter to measure the position quadrature 
X. 

The evolution equation for the total density operator D for the vibrational degree of freedom and the internal states 
is determined by Hamiltonian (|]) and by the terms describing the spontaneous emission from the level |+) responsible 
for the fluorescence, 

C-spontD = - (2<7-Dg + - <7 + <7-D - Do+oJ) , (5) 

where k is the spontaneous emission rate. Here we have neglected the recoil and the associated heating of the vibra- 
tional motion. This is reasonable in the Lamb-Dicke limit we have assumed from the beginning, since the associated 
heating rate is given by K(kxo) 2 vibrational quanta per second, which is negligible for a sufficiently small Lamb-Dicke 
parameter kxQ. In practical situations, also other heating mechanisms exist, caused by technical imperfections such 
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as the fluctuations of trap parameters due to ambient fluctuating electrical fields in the ion trap case 0, and due 
to laser intensity noise and beam-pointing fluctuations in the case of far-off resonance optical traps (see Ref. || and 
references therein). We assume the presence of this heating due to trap imperfections, and we describe it with the 
following term in the master equation, characterized by a heating rate 7^ (see ]l2|]) 

C h D = y (2aDa) - a)aD - Da) a) + y (2a) Da - aa) D - Daa)) . (6) 

The heating rate 7^ has not to be too large, in order to stay within the assumed Lamb-Dicke regime. The Lamb-Dicke 
condition also implies that the trapped atom has to be initially prepared in a sufficently cold state, i.e., an effective 
thermal state with, say, a mean vibrational number no ~ 10. This can be obtained with a preliminary Doppler cooling 
stage, which is then turned off at t — and replaced by the proposed feedback cooling scheme. We shall see that our 
scheme is able to further cool the trapped atom, close to the ground state, even in the presence of moderate heating 
processes. 

The resulting master equation for the internal and vibrational degrees of freedom is 

D = C h D - I [H,D] + I (2o-Do+ - a+a^D - Da+a^) . (7) 

Let us now see how to realize the continuous position measurement. It has been recently shown that when excited 
by a low intensity laser field, a single trapped atom emits its fluorescent light mainly within a quasi-monochromatic 
elastic peak [JL3|. The fluorescent light spectrum was measured by heterodyne detection. By improving the technique 
it does not seem impractical to get a homodyne detection of the single-ion fluorescent light. In Ref. it was shown 
how one could achieve such a measurement. Thus, by exploiting the resonance fluorescence it could be possible to 
measure the quantity £ v = (cr_e~ iv + <r+e l(p ) through homodyne detection of the field scattered by the atom along 
a certain direction ||. In fact, the detected field may be written in terms of the dipole moment operator for the 
transition |— ) «-> |+) as || 

fiW(t) = ^_(t), (8) 

where r\ is an overall quantum efficiency accounting for the detector efficiency and the fact that only a small fraction 
of the fluorescent light is collected and superimposed with a mode-matched oscillator. 
As a consequence of @), the homodyne photocurrent will be jl5| 

I(t) = 2r}K{-£ v {t)) c + V^eW , (9) 



where the phase tp is related to the local oscillator, which, since we have assumed the resonance condition A = 0, in 
the present case is provided by the same driving field generating the classical standing wave. The subscript c in Eq. (|^) 
denotes the fact that the average is performed on the state conditioned on the results of the previous measurements 
and is a Gaussian white noise Epj . In fact, the continuous monitoring of the electronic mode performed through 
the homodyne measurement, modifies the time evolution of the whole system, and the state conditioned on the result 
of measurement, described by a stochastic conditioned density matrix Z? c , evolves according to the following stochastic 
differential equation (considered in the Ito sense) 

D c = C h D c - % - [H, D c ] + ~ (2<t_L> c ct + - a+a^D c - D c a+a^) 

+ y/VKt(t) (e-^a^D c + e llfi D c a+ - 2{X V ) C D C ) . 

Since we are considering a strong fluorescent transition it is reasonable to assume that the spontaneous emission 
rate k is large, i.e., k>^. This means that the internal two- level system is heavily damped and that it will almost 
always be in its lower state |— ). This allows us to adiabatically eliminate the internal degree of freedom and to 
perform a perturbative calculation in the small parameter x/k, obtaining (see also Ref. GJ|) the following expansion 
for the total conditioned density matrix D c 

D c = Pc ® \-)(-\-i~(Xp c ®\+)(-\-p c ® \-)(+\X) , (10) 

where p c = Tr e ; D c is the reduced conditioned density matrix for the vibrational motion. In the adiabatic regime, 
the internal dynamics instantaneously follows the vibrational one and therefore one gets information on the position 
dynamics X by observing the quantity E v . The relationship between the conditioned mean values follows from Eq. (nQ) 
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(E v (t)) e = ±{X(t)) c 8m<p. 



(11) 



Moreover, if we adopt the perturbative expression ( |l0| ) for D c in ( |10|) and perform the trace over the internal mode, 
we get an equation for the reduced density matrix p c conditioned to the result of the measurement of the observable 
(S v (t)) c , and therefore (X(t)) c 

x 2 

p c = C h p c - iv \aJa, p c ] - — [X, [X, p c ]] 

Ik 

+ VvX 2 /^m (ie^p c X - %e-**Xp c + 2s\n V {X{t)) cPc ) . (12) 

This equation describes the stochastic evolution of the vibrational state of the trapped atom conditioned to the result 
of the continuous homodyne measurement of the resonance fluorescence. The double commutator with X is typical 
of quantum non-demolition (QND) measurements of the position. However this indirect measurement is not properly 
QND because of the presence of the vibrational bare Hamiltonian Tiva) a mixing the position quadrature with the 
momentum. 



B. The Off-Resonant Interaction 



One can realize an effective indirect measurement of the atomic position also in the opposite limit of large detuning 
between the standing wave field and the internal transition. In fact, when the internal detuning A is very large 
(A ^> v, 5, k, e) the excited level |+) can be adiabatically eliminated: the state of the whole system (atom+standing 
wave mode) can be written as "0+1+) + ^-| — )■> an d adding a constant term KA/2 to the Hamiltonian (Q), the 
corresponding Schrodinger equations will be 

iip + = (A + va^a + 5tfb)ip + + ie sin(fcr + <f>)top- (13) 
= {va) a + Sb^b)^- — iesxn(kx + (f>)b^ip + . (14) 



In the adiabatic limit of very large A we can neglect the time derivative in (13) and put 



ip + ~ — i— sin(fcx + (f))bip- . (15) 

Inserting this equation into (|l4|), one gets an equation for which is equivalent to have the following effective 
Hamiltonian for the vibrational motion of the atom and the standing wave mode alone 

f 2 

H = h6tfb + hi/a t a-h— b^bsm 2 {kx + <j)) . (16) 
If we now set the spatial phase 4> — 7r /4, we can rewrite ( |l6|) as 

H = h(s- ^\ b f b + hvafa - ft^&Wn (2kx) . (17) 

It is clear that in this case it is convenient to choose the frequency u> of the rotating frame so that 6 — e 2 /2A. The 
Hamiltonian ( |l7| ) assumes the desired form when the Lamb-Dicke regime is again assumed so to approximate the sine 
term with its argument, and when the case of an intense standing wave is considered. However, in this case we shall 
not neglect the quantum fluctuations of the standing wave field, and we shall make the replacement b — > [3 + 6, where 
/3 > 1 describes the classical coherent steady state of the radiation mode and b is now the annihilation operator 
describing the quantum fluctuations. One gets 

,2 

H = hva ! a-n— kx(\/3\ 2 + (3*b + (3b^) . (18) 

Shifting the origin along the x direction by the quantity h€ 2 k\(3\ 2 / /S.mv 2 : one finally gets an effective Hamiltonian 
analogous to that of the resonant case (Q) 

H = Tiva)a + HxYX , (19) 
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where now \ = — 4|/3|fcxoe 2 /A and the atomic polarization a x is replaced by the standing wave field quadrature 
Y = (be^ 1 ^ 13 + we l ^ )/2, where (f>p is the phase of the classical amplitude j3. This means that in this case the "meter" 
is represented by the cavity mode, and that an effective continuous measurement of the position of the trapped atom 
is provided by the homodyne measurement of the light outgoing from the cavity. This measurement allows in fact to 
obtain the quantity Y v — {aeT %v +a' f e lip )/2, which is analogous to the quantity E v of the previous Section. Therefore, 
all the steps leading to Eq. ( |l2| ) in the preceding subsection, can be repeated here, with the appropriate changes. 
In this non-resonant case, D now refers to the density matrix of the system composed by vibrational mode and the 
standing wave mode and the spontaneous emission term in the master equation (Q) has to be replaced by the formally 
analogous term describing damping of the standing wave mode due to photon leakage. This is equivalent to interpret 
the parameter k as a cavity mode decay rate in this case and to replace ct_ with b, <r+ with b', and X v with Y v 
in Eqs. (]?]), (||), @ and (pX|). It is again reasonable to assume that the standing wave mode is highly damped, i.e. 
k x, so that it is possible to eliminate it adiabatically. The perturbative expansion ( |l0| ) now becomes 

D c = (p c -^Xp c x) ®|0)(Q|-^(Xp c O|l)(Q|-p c X®|0)(l|) 

\ K z J K 

2 2 

+ X-X Pe X®\l){l\--£-=(X 2 p D ® \2)(Q\+ Pc X 2 ® |0)(2|) , (20) 

K K Z V2 

where \n), n = 0,1,2, are the lowest standing wave mode Fock states. Using this adiabatic expansion and tracing 
over the standing wave mode, one finally gets exactly Eq. (|l2]), describing the reduced dynamics of the vibrational 
mode conditioned to the result (X(t)) c of the continuous position measurement. 

III. THE FEEDBACK LOOP 

We are now able to use the continuous record of the atom position to control its motion through the application of 



a feedback loop. We shall use the continous feedback theory proposed by Wiseman and Milburn |17 |. 

One has to take part of the stochastic output homodyne photocurrent I(t), obtained from the continuous monitoring 
of the meter mode, and feed it back to the vibrational dynamics (for example as a driving term) in order to modify 
the evolution of the mode a. To be more specific, the presence of feedback modifies the evolution of the conditioned 
state p c (t). It is reasonable to assume that the feedback effect can be described by an additional term in the master 
equation, linear in the photocurrent I(t), i.e. fl7|| 

[Pc(«)] / 6=^^^c(«) ) (21) 

1 vx 

where r is the time delay in the feedback loop and K. is a Liouville superoperator describing the way in which the 
feedback signal acts on the system of interest. 

The feedback term (^l|) has to be considered in the Stratonovich sense, since Eq. (^l|) is introduced as limit of a 
real process, then it should be transformed in the Ito sense and added to the evolution equation (|l2|). A successive 
average over the white noise £(t) yields the master equation for the vibrational density matrix p in the presence of 
feedback. Only in the limiting case of a feedback delay time much shorter than the characteristic time of the a mode, 
it is possible to obtain a Markovian equation p7| , p^| , which is given by 

p = C hP - iv [o+o, p] - \- [X, [X, p\] + K (ie^pX - ie^Xp) + tt^TT.P- ( 22 ) 

The second term of the right hand side of Eq. ( |22| ) is the usual double-commutator term associated to the measurement 
of X; the third term is the feedback term itself and the fourth term is a diffusion- like term, which is an unavoidable 
consequence of the noise introduced by the feedback itself. 

Then, since the Liouville superoperator K, can only be of Hamiltonian form []l7f , we choose it as K,p = —ig [P, p] /2 
||l6|| , where P = (a — a^)/2j is the adimensional momentum operator of the trapped particle and g is the feedback 
gain related to the practical way of realizing the loop. One could have chosen to feed the system with a generic 
phase-dependent quadrature; however, it is possible to see that the above choice gives the best and simplest result 
p6| . Using the above expressions in Eq. ( |22| ) and rearranging the terms in an appropriate way, we finally get the 
following master equation: 
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where 



p = ~(N + 1) (2apa) — a) ap — pa) a) + — N (2a) pa — aa) p — paa) 



—M (2a) pa) - a) 2 p - pa) 2 ) - —M* (2apa - a 2 p - pa 2 
— sin tp [a 2 — a) 2 , p\ — iv [a) a, p] , 



(23) 



r = — gsmtp : 
N = -_ * 



M 



<7 sm tp 
1 



ir]x 2 /K 



1 

2 : 



<7 sm 93 



XT 
4k 



Arjx 2 /K 



— cot (fl . 
2 ^ 



(24) 
(25) 

(26) 



Eq. (J2^) is very instructive because it clearly shows the effects of the feedback loop on the vibrational mode a. The 
proposed feedback mechanism, indeed, not only introduces a parametric driving term proportional to g sin tp, but it 
also simulates the presence of a squeezed bath pj[ , characterized by an effective damping constant T and by the 
coefficients M and N, which are given in terms of the feedback parameters p6[ . An interesting aspect of the effective 
bath described by the first four terms in the right hand side of (|2^) is that it is characterized by phase-sensitive 
fluctuations, depending upon the experimentally adjustable phase tp. This master equation preserves the positivity 
of p provided that the condition \M\ 2 < N(N + 1) is satisfied p9[, as it can be checked in the present case. In 



fact, under this condition it is always possible to find a unitary transformation transforming Eq. 
equation manifestly of the Lindblad form. 



Q23|) into a master 



IV. THE STATIONARY SOLUTION 



Because of its linearity, the solution of Eq. (|2^) can be easily obtained by using the normally ordered characteristic 
function [fl9| C(X,X*,t). The partial differential equation corresponding to Eq. (p3|) is 



dt + + ii^j Xd x + (I - x *dx* + I sin tp (\d\* + A : <\ 1 } C( \. A : . / > 
TW|A| 2 + ^M+|sin^ (A*) 2 + f^M* + ^ S mtp) \ 2 } C(\. X' . I) . 



(27) 



The stationary state is reached only if the parameters satisfy the stability condition that all the eigenvalues have 
positive real part, which in the present case is achieved when gsimp < 0. In this case the stationary solution has the 
following form 



where 



C(A, A*, 00) = exp 



-C|A| 2 + i M (A*) 2 + ^*A 2 



N(g 2 sin 2 tp + Av 2 ) + g sin tp (2i/Im{M} - g sin tpRe{M}) + g 2 sin 2 tp/2 



Av 2 



p = T 



(N + 1/2)3 sin tp + TRe{M} + 2vIm{M} 
4^2 



ii^P [Re{M }-(N + 1/2)] 



(28) 



(29) 



(30) 



This means that in general, the stationary state is a generalized gaussian state. However, in practical situations, 
the stationary state assumes a much simpler form. In fact, the vibrational frequency v is usually much larger than 
the heating rate 7/j and the feedback parameter g, and in this limit, Eqs. (29) and ( pfj| ) become ( w N and /1 w 
respectively. Using Eq. (Ka), one has 



C(A,A*,cx)) = exp [— 7V|A| 2 ] 



(31) 



G 



that is, the stationary state is an effective thermal state with mean vibrational number N. This can also be seen from 
the fact that in the large v limit, one can consider the master equation ( p3| ) in the frame rotating at the frequency 
v, and neglecting the rapidly oscillating terms, one ends up with a thermal master equation given by the first line of 
Eq. (p3|), whose steady state is just the thermal state with mean phonon number N . 

The expression for N given by Eq. ( pjj ) shows that it is convenient to choose <p — —tt/2 to get the smallest possible 
values for N. In this way the stability condition is also automatically satisfied. Then, the minimum value for the 
stationary mean vibrational number N can be obtained by minimizing it with respect to the feedback gain g: the 

1/2 

optimal value for g is given by g = 4 ((7/1 + x 2 /4k) i]x 2 /4k) and the corresponding minimum value of ./V is 



This expression shows the best cooling result achievable with the present feedback scheme. One has that when the 
heating rate is negligible with respect to the "QND coupling" parameter % 2 /4k, the final vibrational number ./V is 
limited only by the efficiency of the homodyne measurement 77. In particular this fact is true in the optimal case 
where all the technical sources of heating are eliminated which means jh = in all the above expressions. Therefore, 
the present scheme is able to achieve cooling to the motional ground state especially in the off-resonant scheme, in 
which the particle position measurement is realized through homodyning the radiation exiting the cavity. In fact, 
in this case the homodyne efficiency can be very close to one. This is shown in Fig. 1, where we have sketched the 
phase space uncertainty contours obtained by cutting the Wigner function corresponding to Eq.(|2S|) at \j \fe times its 
maximum height. We see that the feedback produces a relevant contraction of the uncertainty region, which becomes 
almost indistinguishable from the region corresponding to the motional ground state (inner dotted line in Fig. 1). The 
outer dashed line corresponds instead to the initial thermal state with mean vibrational number n = 10, prepared 
by the Doppler pre-cooling stage. 

In the resonant case in which the effective position measurement is realized through the homodyne measurement of 
the fluorescence, the measurement efficiency is much lower and ground state cooling becomes very difficult to achieve. 
However, this position measurement scheme based on fluorescence becomes necessary when one cannot extract the 
light out of a cavity, such as in Ref. [Q~3|, and one has to use the fluorescent light. 



We have proposed a feedback scheme able to achieve significant cooling of the motional degree of freedom of a 
trapped particle. The method is based on an effective continuous measurement of the particle position which can 
be realized in two different ways: by homodyning either the fluorescence of a strong transition or directly the light 
exiting the cavity which is coupled to the trapped particle. When the efficiency of the homodyne detection is close 
to one, the method is able to achieve ground state cooling. It is interesting to note that, even if only the particle 
position is measured and the feedback is chosen in order to decrease position fluctuations, the scheme provides a 
phase-space localisation for all quadratures. This is essentially due to the fact that the bare atom Hamiltonian hva^a 
mixes the dynamics of the atomic position and momentum, so that the continuous homodyne measurement actually 
gives informations on both quadratures. This model shares some peculiarities with that one we have proposed in pc(| 
to cool the vibrational motion of a macroscopic mirror of an optical cavity. The present application to a trapped 
atom stresses the versatility of the methods using feedback loops in systems characterized by the radiation pressure 
force in controlling thermal noise. 

It is also possible to see that the proposed scheme is not able to reduce the noise below the quantum limit, i.e. 
the stationary state variance of a generic motional quadrature cannot be made smaller than 1/4. This is again a 
consequence of the free vibrational Hamiltonian hva^ a (it is in fact possible to get position squeezing in the limiting 
case v = pl[ ) . Actually, position squeezing can be obtained by considering a suitable modification of the present 



As concerns the specific way in which a particular feedback Hamiltonian could be implemented, the important point 
is to be able to realize a term in the feedback Hamiltonian proportional to momentum. This is not straightforward, 
but could be realized by using the feedback current to vary an external potential applied to the atom without altering 
the trapping potential. On the other hand, shifts in the position (being strictly equivalent to a linear momentum 
term in the Hamiltonian), are achieved simply by shifting all the position dependent terms in the Hamiltonian, in 
particular the trapping potential. Alternatively, the use of laser pulses could be useful as well, since, using a typical 
laser cooling scheme, the light can exert on the atom a force proportional to its momentum. 




(32) 



V. CONCLUSIONS 



scheme p2[ . 
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As we have already remarked, in principle the model could be extended to the three dimensional case. As concerns 



the model discussed in Sec. [I A, one should consider three different internal transitions, each one coupled with a 



vibrational degree of freedom, resonant with three orthogonal standing waves. For the off-resonant case presented in 



Sec. [IB, one should only consider three orthogonal standing waves far from resonant transitions. 

In conclusion, although the implementation of the presented cooling method via feedback could be an hard task, 
it can be useful whenever the use of resolved sideband cooling is impractical. The possibility of having an alternative 
way to cool trapped particles is particularly interesting for quantum information processing applications, because it 
may happen that the requirement of having two highly stable internal states for quantum logic operations and a good 
internal transition for sideband cooling cannot be simultaneously satisfied. With this respect other cooling strategies 
have been recently proposed, as for example the use of sympathetic cooling between two different species of ions P,p3[. 
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FIG. 1. Phase space uncertainty contours obtained by cutting the Wigner function of the stationary state at l/\/e times 
its maximum height. The dashed line refers to the initial thermal state with mean vibrational number no = 10; the solid line 
refers to the steady state in the presence of feedback with \ — 4 kHz, g = 0.375 kHz, v = 1 MHz, 7>, = 10 Hz, k = 40 kHz, 
t/ — 0.9, tp = — 7r/2. Notice that feedback provides cooling very close to the ground state (the corresponding contour is given 
by the dotted line) . 
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